Three-Dimensional Analysis of Wakefields Generated by 
Flat Electron Beams in Planar Dielectric-Loaded Structures 
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An electron bunch passing through dielectric-lined waveguide generates Cerenkov radiation that 
can result in high-peak axial electric field suitable for acceleration of a subsequent bunch. Axial 
field beyond Gigavolt-per-meter are attainable in structures with sub-mm sizes depending on the 
achievement of suitable electron bunch parameters. A promising configuration consists of using pla- 
nar dielectric structure driven by flat electron bunches. In this paper we present a three-dimensional 
analysis of wakefields produced by fiat beams in planar dielectric structures thereby extending the 
work of Reference [A. Tremaine, J. Rosenzweig, and P. Schoessow, Phys. Rev. E 56, No. 6, 7204 
(1997)] on the topic. We especially provide closed-form expressions for the normal frequencies and 
field amplitudes of the excited modes and benchmark these analytical results with finite-difference 
time-domain particle-in-cell numerical simulations. Finally, we implement a semi-analytical algo- 
rithm into a popular particle tracking program thereby enabling start-to-end high-fidelity modeling 
of linear accelerators based on dielectric-lined planar waveguides. 

PACS numbers: 29.27.-a,41.20. Jb,41.60.Bq 



I. INTRODUCTION 

Next generation multi-TeV high-energy-physics lep- 
ton accelerators are likely to be based on non- 
conventional acceleration techniques given the limitation 
of radio-frequency (rf) normal-conducting [H and super- 
conducting 2] structures. Non-conventional approaches 
based on the laser plasma wakefield accelerator have re- 
cently demonstrated average energy gradients one order 
of magnitude higher than those_ possible with state-of- 
the art conventional structures ^3|- The applicability of 
laser-driven techniques to high-energy accelerators is cur- 
rently limited as attaining luminosity values similar to 
those desired at the International Linear Collider would 
demand a laser with power approximately four orders of 
magnitude larger than the most powerful lasers currently 
available [3]. Another class of non-conventional accel- 
erating techniques includes beam-driven methods which 
rely on using wakefields produced by high charge drive 
bunches traversing a high-impedance structure to accel- 
erate subsequent witness bunches Such an approach 
has the advantage of circumventing the use of an external 
power source and can therefore operate at mm and sub- 
mm wavelengths. Structures capable of supporting wake- 
field generation include plasmas [HI , and dielectric- loaded 
waveguides [l3|. The possible use of plasma-wakefield 
accelerators (PWFAs) as the backbone of a a multi-TeV 
electron-positron linear collider is limited by plasma ions 
motion due to the intense electromagnetic field of the 
bunch . Dielectric wakefield accelerators (DWFAs) are 
not prone to similar limitations. 

In this paper we concentrate on the colincar DWFA. In 
such a configuration a highly-charged drive bunch prop- 
agates through a dielectric-lined waveguide (DLW) and 



excites an electromagnetic wake [i-li^. A delayed wit- 
ness bunch moving on the same path as the drive bunch 
can experience an accelerating field. 

Recent experiments [13] confirm that DLW can sup- 
port accelerating fields in excess of a GV/m thereby mak- 
ing DWFA a plausible candidate for the next generation 
high-energy-physics linear accelerators [ij] or compact 
short- wavelength free-electron lasers {l^ . 

In cylindrically-symmetric DLW structures, the elec- 
tric field amplitude of the wakefield is approximately in- 
verse proportional to the aperture radius. Given the lin- 
ear charge-scaling of the field, peak electric field ampli- 
tude of the order of Gigavolt-per-Meter can be obtained 
by either using high-charge drive beams (> 100 nC) in 
mm-sized DLW structures or by focusing sub-nC bunches 
in micron-sized DLW structures d, [111, jl^ • 

To date, cylindrically-symmetric DLW's have been ex- 
tensively studied both theoretically ^9,, 16- 18] and ex- 
perimentally [H-dl] . Since the angular divergence of the 
beam also increases during focusing it is hard to maintain 
a low transverse size of the beam over a long propaga- 
tion distance. Therefore, the design of the DLW struc- 
tures must compromise between a small transverse size 
to maximize the intensity of the wakefield, and a longer 
interaction length to maximize the energy gain of the test 
charge. An appealing solution consists of using flat elec- 
tron beams (cTx ^ (Jy) passing through slab-geometry 
DLWs [m, [23i - |25j . This possibility has become more at- 
tractive since the recent advances toward generating fiat 
beams directly in photoinjectors [2^ [27j. 

In the following Sections we present an analysis of the 
generation of wakefields in rectangular DLW structures 
excited by drive beams with arbitrary three-dimensional 
charge distributions. The main goal of this paper is 
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to extend the formalism introduced in Ref. [llj by in- 
cluding all types of modes excited in a planar DLW. 
This paper is pedagogical in the sense that it is about 
the method of derivation rather than the results, which 
have been derived in previous papers, most extensively 
in Ref. [IS]- However, this paper provides close- form for- 
mulae for the eigenfrequencies and electromagnetic fields 
excited in planar DLWs. These analytical results are 
benchmarked against three-dimensional finite-difference 
time-domain (FDTD) simulations. Finally, the model 
is implemented in a popular particle-in-cell (PIC) beam 
dynamics program. The latter provides a fast and high- 
fidelity model enabling start-to-end simulation of DLW- 
based linear accelerators. 



II. WAKEFIELD GENERATION 

The geometry of the problem analyzed in this paper is 
depicted in Fig. [TJ Transient effects resulting from the 
injection of the electron bunch in the structure are not 
included (the structure is assumed to be infinitely long) 
and the drive bunch is taken to be ultra-relativistic with 
its Lorentz factor 7 much greater than unity. The build- 
ing block of a real drive bunch is assumed to consist of a 
linear charge distribution oriented along the x-axis which 
moves in the z-direction with velocity v and has offset ?/o 
in the vertical direction. To simplify, it is also assumed 
that the linear charge distribution is symmetric in the 
x-coordinate and it vanishes at the ends of the dielectric 
structure. Under these assumptions the charge distribu- 
tion can be written as a Fourier series in the horizontal 
direction 11^ 

p[x, y, z) = ^ A„ cos{kx^rnx)5[y - yo)5[z - vt), (1) 

where Am is a constant and each term is indexed by 
the integer m — 0, 1,--- defined such that k^.m = 
(2m + 1)^. In the charge- free vacuum region the elec- 
tromagnetic Wakefield satisfies the wave equation: 



V^-l— |E = 



(2) 



In this paper we consider only the propagating modes 
with harmonic longitudinal and temporal dep endencies 
of the form E(a;, y, z, t) = E(a;, y)e*("*-'=-^) [lit HI where 
kz ~ k.z. In addition to the propagating- mode solutions, 
Eq. [2]when subject to the boundary conditions also ad- 
mits evanescent modes with imaginary eigenfrequencies 
w [2^. Since these evanescent modes are present over 
short distances behind the drive bunch, they do not con- 
tribute to the beam dynamics of a subsequent witness 
bunch. These short-range fields are consequently ignored 
in the remaining of this paper. 

In the ultra-relativistic regime the synchronism con- 
dition is achieved because the wakefield phase veloc- 
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FIG. 1: Overview of the DLW structure geometry (left) and 
transverse cross section (right). A rectangular-shaped drive 
beam is displayed in blue. 



of the witness beam. 



Since the transverse wave vector 
k/j — the wave equation [2] can 



k± = ^B-kl 
be written only in terms of the partial derivatives with 
respect to the transverse coordinates: 



ViE(a;,y) =0. 



(3) 



The symmetry of the drive-beam charge distribution and 
the boundary conditions at a; = ±Lx/2^ determine (up to 
constants) the analytical expression of the fields. In the 
vacuum region the electric field components are simple 
combinations of trigonometric and hyperbolic functions: 



Ex (X sv[\{kx^r,ix)cos\v{kyy), 
Ey cx cos{kx^„ix)smh{kyy), 
Ez oc cos{kx^mx) cos\i(kyy). 



(4) 



The axial field Ez associated with this set of solution 
is symmetric with respect to the horizontal axis [i.e. 
Ez{x, —y, z) = Ez{x, y, z)] we henceforth refer this set to 
as "monopole" modes. Similarly, a set of fields with an 
antisymmetric axial field [Ez(x, —y, z) = —Ez{x,y,z)] is 
obtained by substitution of sinh(fcyy) with cosh{kyy) and 
vice versa. This latter set is termed as "dipole" modes 
in the remaining of this paper. 

Inside the dielectric, the transverse wave number can- 
not be neglected 



k I — 



(er - 1) > 0, 



(5) 



where is the relative electric permittivity of the dielec- 
tric medium. The boundary conditions at y = ±6 and 
at a: = ±Lj^/2 determine the trigonometric form of the 
field expressions. Inside the dielectric region, there is no 
distinction between monopole and dipole modes and the 
electric field components are given by 



Ex cx sm{kx^m,x)cos[ky{b ~ y)], 
Ey cx cos{kx^mx)sm[ky{b ~ y)], 
Ez cx cos{kx^mx)cos[ky{b - y)]. 



(6) 



ity 



c is very close to the velocity 



The corresponding expressions for the magnetic field in 
the vacuum and dielectric regions can be easily obtained 
following a similar prescription. 



3 



Inspection of the z-component of the fields indicates 
that the normal modes cannot be categorized in the usual 
transverse electric or transverse magnetic sets. The rea- 
son is that the separation surface between the dielectric 
and vacuum regions is in the x — z plane unlike the case 
of the uniformly-filled waveguides where this surface is in 
the transverse x — y plane. 

Therefore, it is natural to categorize the modes depend- 
ing on the orientation of the fields with respect to the 
dielectric surface. Following the definition introduced 
in Ref. [l^ we classify the mode as Longitudinal Sec- 
tion Magnetic (LSM) and Longitudinal Section Electric 
(LSE) modes corresponding respectively to the case when 
the magnetic and electric field component perpendicular 
to the dielectric surface vanishes. In the configuration 
shown in Fig. [1] LSE and LSM modes correspond re- 



spectively to Ey 
interface. 



and = at the vacuum-dielectric 



III. DISPERSION EQUATIONS & 
ELECTROMAGNETIC WAKEFIELDS 

A. Longitudinal Section Magnetic (LSM) modes 

To obtain the normal mode frequencies it is convenient 
to use the hertzian potential vector method [1^. In the 
source-free region, magnetic and electric fields can be ex- 
pressed in terms of the vector potential 11: 



dxdy 



E, 



^'^^ dy 



Hy=0 



Hz — iekzC 



(9) 



dx 



The expression for ijjf. that satisfies Eqs.|ni[Sl and [7] is 
given by: 



■0e 



Acos{kj:^„ix) sinh(fca;,my), < y < a 
B cos(fcj. „a;) cos [ky{h — y)] ^ a < y < b 



(10) 



where A and B are constants. The boundary condi- 
tions completely determine the eigenfrequencies of the 
LSM modes. Since (cx eipe) and E^ (oc ^^) are con- 
tinuous at y — a the constants A and B can be eliminated 
giving the dispersion equation 



coth(fc^^,„a) cot [ky{b — a)] = 



^rkx 



(11) 



Therefore, for each discrete value of kx^m there is an 
infinite set of discrete fcj,_„ values where n is an integer. 
So, the eigenfrequencies are indexed by the integer couple 
(m, n) and verify 



H = iLoeV X n, 



E 



/c^n -h v(v • n). 



(7) 

(8) 



where 11 is the solution of the homogeneous Helmholtz 
equation in the vacuum region. To determine the fre- 
quencies and the fields for the modes when Hy = (LSM) 
it is convenient to choose 11 = He = '!/'e(a;, y)e*('^*^''^^)y 
where the the scalar function 4^eix, y) must be a solution 
of (V5_ + k\)Tpe = 0. From Eqs.[7]and[5]all fields can be 
expressed in terms of the unknown function ipe^x^y): 



k^ - 



(12) 



It is worthwhile noting that the boundary conditions do 
not completely determine either the function or the 
fields. The reason being that, up to this point, the field 
source terms (p, j) were not taken into account although 
their symmetries were invoked. Still, it is straightforward 
to show that all fields (and also 'ipe) associated to a certain 
mode depend on a common normalization constant, the 
amplitude -Eo;m,n, which remains to be determined. The 
expressions of the fields are given by 
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Ex 



»go;m,nfc»,m sm{kx^rnx) cosh(fcj:,„y) < y < a 

iEo.,^:nk.,^ cosh(fc..^_a) 3in(fc,_^a;) sin [fc,,„(6 - y)] a<y<6 



sin[fc^_„(6— a)] 



fefSr (^^,"1 + fcz) cos(fcj;,„ia;) sinh(A;^,„y) < y < a 

^ CnCsiniCnffe-S (fe^ + ^x,m) cos(fc^,„a;) cos [fcj,,„(6 - y)] a<y<b 



Eo-rn,n COs{kx,mX) COSh{kx,my) <y < a 

^Eoim,n ^Z[ky!7(b'-a)] '^°^i^^,rnx) sin [ky,n{b ~ v)] a<y<b 



H. 



iEo;rn,r.k^ec ^Q^fJ^^^^^^^ sinh(fca;^„y) 



< 2/ < a 



^,m,n ■y_ ^Eo■,^^^k.ec 2°^^"^^^-^^^] COsik,,^x) COS [fc^,„(& - y)] a<y<b 



I fc^.n sin 



< y < a 

'[ky'^n(b-a)] sin(fecc,ma:) COS [fcj^,n(6 - y)] a<y<b 



(13) 



B. Longitudinal Section Electric (LSE) modes 

The case of the LSE modes {Ey = 0) can be treated 
in the same way as the LSM modes. The hertzian vector 
electric potential is replaced by a vector magnetic poten- 
tial Uh which is related to the fields: 



E = —iujfiV X Ilh 

H = /c^Hh + V(V • Hh) 



(14) 
(15) 



As in the previous case, it is convenient to factor out 



the t and z-dependencies of the Ilh and to define a 
scalar function iph{x,y)'- Hh = iph{x,y)e''^"*~'^''^^y. It 
is straightforward to derive the dispersion equation for 
the LSE modes 



coth(A;a;a) cot [ky{b — a)] = — t^, (16) 

ky 



and the expressions for the electromagnetic-field compo- 
nents: 
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iEo 



<y < a 



sinTfcl^llfe-a)] ^Mkx,mx) sin - y)] a<y<b 



K — f) 



Eo COs{kx^mX) COSh{kx.rny) 



En 



'o;m,« s'in [ fc^'^rb-a) ] cos(fca;,ma;) sin [fcy,„(6 - ?/)] a <y <b 



< y < a 



Her 



fcj^^e"''" cos(fc2:,„a;) s\nh{kx^r,iy) <y < a 

'^°7:fc'" ZtylT-a)] CQs(fcx,ma:) cos [fc^,„(6 - 7/)] a<y<6 



y,m,n 



sin(fc^,ma;) cosh{kx^ray) 



'''^"'"^d^lkZc^''^ sinltj.lira)] sin(fcx,ma;) COS [fcj,,„(6 - y)] a<y<b 



■ siTi{kx^mx) sinh(fc; 

,^nky cosh(A: 



y) <y <a 
I 



(17) 



IV. WAKEFIELD AMPLITUDES 

Although we have obtained the general expression for 
the electromagnetic field, the constants i?o;m,n still re- 
mained to be determined. An often used method to find 
the constants Eo-m.n in Eqns. [13] and 1171 is to evaluate 
the Green function of the wave equation with sources in- 
cluded. The latter is usually a cumbersome procedure 
and we instead chose to determine i?o;m,Ti's based on en- 
ergy balance considerations. The electromagnetic energy 
stored in the DLW equals the mechanical work performed 
on the drive bunch. For a linear wakefield, the decelerat- 
ing field Ed acting on a drive point-charge is half of the 
wakefield amplitude, i.e. Ed = as a consequence of 

the fundamental wakefield theorem f39 |. . Suppose the 



drive charge moves over an infinitely short distance Sz. 
The work performed by the wakefield on the drive charge 
should equal the energy stored in the field 



Ez;m,niz = Vt) 



dxdySz 



Ue'^ + i-i^H^) dxdySz, 



(18) 



where the integration extends over the transverse plane. 
It is straightforward to evaluate the full expressions of 
the wakefield amplitudes in terms of the wave numbers 
k^ „ and fc„,„ from Eqns. [T51 1131 and 1171 



E, 



LSM 
0;m,n 



1 



2eQ sinh(2fcx 



\ra COsh{kx^rnyo) 



,_al tr- cosh^(fcx.„a) r 6_Q /-, „ \ _ sin[2fcy,„ (b-g)] /-, _ \ 



E, 



LSE 



\m COSh{kx^rnyo) 



2fex,™ sin^ [fey, „(b-a)] \ 2 \^ J 4ky,r, [j^ ) } 



r 



(19) 



The field amplitudes for the dipole modes can be smh{kx.m'U') for cosh(fca;^mu) where u takes on the val- 
obtained from the previous equations by substituting 
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ues of a and yo. As expected, the wakefield amplitude 
scales linearly with the drive bunch charge and inverse 
proportionally with the transverse size of the structure. 
To this point this model is fully three-dimensional and 
the only limitation stems from the assumed symmetry of 
the transverse drive charge distribution with respect to 
the y axis. 
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FIG. 2: (color online) Longitudinal (a) and transverse sections 
(b) of the DLW model implemented in vorpal. The rectan- 
gular box delimited by green blocks represents the 3D com- 
putational domain used in the simulations. Grey, green and 
purple blocks respectively stand for perfectly-matched layer 
(PML), perfectly-conducting boundary (PCB), and dielectric 
material (with associated relative dielectric permittivity e^)- 
The cyan rectangles represent the electron bunch distribution. 
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FIG. 3: Eigenfrequencies [fm,n = Wm,n/(27r)] and associated 
amplitudes for the mode induced in the DLW with param- 
eter listed in Tab. H] The (m, n) LSM and LSE modes are 
respectively shown as red and blue bars . 
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V. COMPARISON WITH 
THREE-DIMENSIONAL FDTD SIMULATIONS 

To evaluate the wakefield associated to a drive bunch, 
an integration over the full three-dimensional continuous 
charge distribution must be performed. In practice the 
integration is replaced by numerical summations of dis- 
crete charge distributions similar to those described by 
Eq. [TJ This process is similar to the charge discretization 
procedure used in standard PIC algorithms |33|. An im- 
portant feature of this model is that the integration over 
x-direction is already performed through the Fourier ex- 
pansion of the charge distribution. For the charge distri- 
butions considered in this paper, only a few Fourier terms 
(< 10) are sufficient to obtain an accurate representation 
of the distribution along the x direction. This is signifi- 
cantly less than the number of grid points in z-direction 
needed by most PIC codes to evaluate the 3D-collective 
effects and external fields. 

The integration over the vertical direction is straight- 
forward and in most cases, when the drive charge distri- 
bution is also symmetric with respect to the horizontal 
axis, the contribution of the dipole modes cancels out. 

Since the phase velocity of the wakefield is the same 
as the velocity of the drive beam, causality principle 
requires that the wakefield vanishes ahead of the drive 
charge. Therefore, the integration over the longitudinal 
direction extends only from the observation point to the 
actual drive charge position. The wakefield assumes the 
form 



E E 

m=0,l,--- n=0,l,- 



p{z')W^,n{z - z')dz' , (20) 



where W stands for any of the electromagnetic field 
components and Wm,n are the corresponding field com- 
ponent associated to the LSMm.„ and LSEm^„ modes 
given respectively by Eq. 1131 and Eq. 1171 The summation 
is performed over all excited modes. 

The results obtained from Eq. [20] are benchmarked 
against simulations performed with vorpal a conformal 
FDTD (CFDTD) PIC electromagnetic solver vor- 
pal is a parallel, object-oriented framework for three di- 
mensional relativistic electrostatic and electromagnetic 
plasma simulation. The DLW model implemented in 
VORPAL is fully three dimensional; see Fig. [51 The 
model consists of the rectangular DLW surrounded by 
perfectly-conducting boundaries (PCBs), The lower and 
upper z planes are terminated by perfectly matched lay- 
ers (PMLs) that significantly suppresses artificial reflec- 
tions of incident radiation [31| . A particle source located 
on the surface of the lower z plane produced macropar- 
ticles uniformly distributed in the transverse plane (x, y) 
and following a Gaussian longitudinal distribution with 
total charge Q described by 



^ix,y,C) 



271(7 zW^Wy 



H 



y,H(f-\y\)H{3a 



ICI), (21) 



where = z — vt, {wy) is the full width transverse 
horizontal (vertical) beam size, and cr^ the longitudinal 
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root-mean-square (rms) length. The function H{...) 
is the Heaviside function. The longitudinal Gaussian 
distribution is truncated at iStTz. Finally, a "particle 
sink" at the upper z plane allows macroparticles to exit 
the computational domain without being scattered or 
creating other source of radiation. 
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FIG. 4: Snapshot of the axial electric field Ez{x = 0,y,z) 
in the mid-plane of a slab DLW for the bunch and structure 
parameters shown in Tab. [H The center of the drive bunch 
is located at « = 27 mm and it is moving in the positive 
z-direction. The field was obtained from vorpal simulations. 

In order to precisely benchmark our model, a DLW 
which supports both LSM and LSE modes is chosen. The 
parameters of the structure and driving bunch are gath- 
ered in Tab. |T] and the frequency and the amplitude as- 
sociated to the first few modes appear in Fig. [3] The 
drive-bunch energy is set to f — 1.0 GeV consistent with 
the ultra-relativistic approximation used in the analytical 
model. 



TABLE I: Parameters of the DLW structure and drive bunch 
used for benchmarking our theoretical model with VORPAL 
simulations. 



parameter 


symbol value unit 


vacuum gap 


a 


2.5 


mm 


height 


b 


5.0 


mm 


width 




10.0 


mm 


relative permittivity 




4.0 




bunch energy 


£ 


1 


GeV 


bunch charge 


Q 


1.0 


nC 


rms bunch length 




1.0 


mm 


bunch full width 


Wx 


6.0 


mm 


bunch full height 


Wy 


4.0 


mm 



The longitudinal component of the electric field sim- 
ulated with VORPAL is shown in Fig. 2] as a two- 
dimensional projection in the y — z plane. 
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FIG. 5: Comparison of the electromagnetic field components 
calculated with our model (blue dashed line) and simulated 
with VORPAL (red solid lines). The fields are computed on a 
line parallel to the z axis with transverse offset a; = 3.0 mm 
and y = 2.0 mm. The center of the drive beam is at 
z = 27.5 mm. The DLW and bunch parameters are the one 
displayed in Tab. HI 



The comparison between theoretical calculation and 
VORPAL simulations are shown in Figs. [5] and [51 The 
fields in Fig. [5] are evaluated as a function of the axial 
coordinate z at a given transverse location (x = 3.0 mm, 
y = 2.0 mm) which corresponds to the upper-left corner 
of the charge distribution. Figure [6] displays the electro- 
magnetic field evaluated as a function of the horizontal 
[Fig. |n]( left)] and vertical [Fig.|n](right)] transverse coor- 
dinate at a given axial location z = 18.6 mm correspond- 
ing to the minimum shown in Fig. [5] (i.e. maximum 
accelerating field). All plots show a decent agreement 
(relative discrepancy < 25 %) between our model and 
VORPAL simulations. The notable disagreement observed 
in Fig. [5] for the transverse fields in the vicinity of the 
driving charge (i.e. z ~ 27.5 mm) is rooted in the absence 
of velocity fields in our model (only radiation field con- 
tributes to the Wakefield) while the VORPAL simulations 
include both velocity and radiation fields. In fact given 
the bunch charge and duration, the amplitude of the ve- 
locity field can be evaluated by convolving the charge 
distribution with the electric field generated by an ultra- 
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FIG. 6: Comparison of the electromagnetic field components 
calculated with our model (blue dashed line) and simulated 
with VORPAL (red solid lines). The fields are computed a a 
given longitudinal position z — 18.6 mm. The left-column 
plots display the dependence on the horizontal coordinate x 
ai y — 0.043 mm while the right-column plots show the de- 
pendence on the vertical coordinate y ai x = 0.043 mm. The 
DLW and bunch parameters are the one displayed in Tab. U 



relativistic particle E(r, ^) — q/ {2n€o){r± /r^)S{() where 
e and eo are respectively the electronic charge and vac- 
uum permittivity, and = (x, y) and C are respectively 
the transverse and longitudinal coordinates of the obser- 
vation point referenced to the particle's location. Such a 
convolution with Eq.[2l]yields the transverse electric field 
components at {x,yX) — (wxjWyjO) to be ~ 1 MV/m 
for the bunch parameter listed in Tab. HI The latter value 
is in agreement with the observed difference between the 



VORPAL and theoretical models; see 
nents in Fig. |6l 



and Ey compo- 



vanishes when — > oo (so that — and to = 0) . In 
such a case the field amplitudes associated to the LSM 
and LSE modes are respectively 



E, 



LSM 
0;0,n 



4A 



a + 



sin-^ \ky_ri {b—a)\ 



and E^.l^, 



0, (22) 



where A is charge per unit length in the horizontal di- 
rection. The LSE modes are suppressed and the latter 
equation is in agreement with the results of Ref. IjJ . 

This limit case can be practically reached by using 
structures with large aspect ratios {Lx 3> Ly) driven 
by flat electron beams (tailored such that Gx ^ <Jy)- 
Flat beams can be produced in photoiniectors by using 
a round-to- fiat beam transformation |33l . Is^ . In such a 
scheme, a beam with large angular-momentum is pro- 
duced in a photoinjector [35l |. Upon removal of the an- 
gular momentum by applying a torque with a set of skew 
quadrupole, the beam has its transverse emittance repar- 
titioned with a tunable transverse emittances ratio [36|. 
Flat beams with transverse sizes of ss 100 ^.m and aspect 
ratio of ~ 20 have been produced, even at a relatively- low 
energy of 15 MeV [13]. 



Vorpal 
-Theory 3D 
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FIG. 7: Longitudinal electric field for a DLW with Lx 3> 



mm and Wy = 



transverse shape with full width ~ 3.0 
150 yLxn, and a Gaussian longitudinal distribution with rms 
length (Tz = 50 /im. The theoretical model based on Eq. [211 
(dashed blue line) is compared with VORPAL simulations (red 
line) . 



VI. THE TWO-DIMENSIONAL LIMIT 

An interesting limiting case occurs when Lx 3> Ly. In 
this "two-dimensional limit" , the dependence of the fields 
on the horizontal coordinate x is weak and completely 



Based on the previous experience in producing fiat 
beams and preliminary simulation of the Advanced Su- 
perconducting Test Accelerator (ASTA) currently in con- 
struction at Fermilab It is reasonable to consider a 
3-nC flat beam generated from the ASTA photoinjector 
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to have parameters tabulated in Tab. [IT] when acceler- 
ated to 1 GeV. Considering a structure with a = 100 ^m, 
b = 300 /im and e — 4.0 would yield a maximum axial 
Wakefield amplitude of 300 MV/m; see Fig. H The 
fundamental LSM frequency is /o,o — 193 GHz and the 
nearest LSE mode amplitude is ~ 136 times lower. In 
Fig. [7] the Wakefield is computed with the asymptotic 
limit provided in Eq. [25] and is in excellent agreement 
with the FDTD simulations. 



TABLE II: Parameter of the DLW structure and drive bunch 
used for benchmarking our theoretical model with vorpal 
simulations. 



this distribution is symmetric the Fourier coefficients A„ 
from Eq. 1 are given by: 



Am - 



2q{xo) COs{kx,rnXo) 



(23) 



In practice the linear charge distribution may not be 
"exactly" symmetric. In this case the charge at an ar- 
bitrary grid point Xq is set at the average of the values 
given by the charge deposition algorithm at positions xq 
and —xq. 



parameter 


symbol 


value 


unit 


vacuum gap 


a 


100 


^m 


height 


b 


300 


^m 


width 




10.0 


mm 


relative permittivity 




4.0 




bunch energy 


£ 


1 


GeV 


bunch charge 


Q 


3.0 


nC 


rms bunch length 




50 


/im 


full (rms) bunch width 


VJx {(Tx) 


3 (0.870) 


mm 


full (rms) bunch height 


Wy (ay) 


150 (43.3) 


/im 



VII. IMPLEMENTATION IN A 
PARTICLE-IN-CELL BEAM DYNAMICS 
PROGRAM 




10 15 20 25 

z (mm) 




10 15 20 25 30 

z (mm) 




z (mm) 
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Although FDTD simulations provide important insight 
that can aid the design and optimization of the DLW 
geometry, their use to optimize a whole linear acceler- 
ator would be time and CPU prohibitive. Therefore it 
is worthwhile to include a semi-analytical version of the 
model developed in this paper in a well-established beam 
dynamics program impact-t [s^. Our main motivation 
toward this choice stems from the availability of a wide 
range of beamline elements models. In addition, impact- 
T takes into account space-charge forces using a three- 
dimensional electrostatic solver. The algorithm consists 
in solving Poisson's equation in the bunch's rest frame 
and Lorentz-boosting the computed electrostatic fields 
in the laboratory frame. 

In a typical particle-tracking PIC program, like 
iMPACT-T, an electron bunch is described by a set of 
"macroparticles" arranged to mimic the bunch phase 
space distribution. Each macroparticle represents a large 
number of electrons (typically lO'^ in our simulations). To 
evaluate the electrostatic fields in the bunch's rest frame 
the macroparticles are deposited on the cells of a three- 
dimensional grid. As a result of this charge deposition al- 
gorithm, the initial charge distribution is approximated 
by a set of point charges, located at the nodes of the 
three-dimensional grid. 

The line charge density corresponding to certain y and 
z-coordinates is p{x) = J2xo lixo)S{x — xq) and provided 




5 10 15 20 25 30 5 10 15 20 25 30 

z (mm) z (mm) 



FIG. 8: Comparison of the electromagnetic field components 
obtained from impact-t [with (green dashed line) and with- 
out (blue dashed line) accounting for space charge effects] with 
the fields simulated with VORPAL (red solid lines). The fields 
are computed on a line parallel to the z axis with transverse 
offset X = 3.0 mm and y = 2.0 mm. The center of the drive 
beam is at z = 27.5 mm. The DLW and bunch parameters 
are the one displayed in Tab. [H 

The general expressions of the wakefields at a given 
position have the following form 



W^ = Yfi{x,y;n) i ^ cos[k^{zo - z)]Ti{zo;n) j , 

n \zo>z / 

(24) 

where i = x,y, z, n is the mode index, /^'s are geometry 
dependent functions of the type described in Eq. [5] and 
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FIG. 9: Longitudinal phase spaces [Sz = z — {z),5l3j = j3^ — £/{mc)] snapshots as a 3-nC electron beam with parameters 
displayed in [Til propagates in a DLW structure at (z) ~ 2.3 cm (a, c) and (z) ~ 22 cm (b, d) from the structure's entrance. 
Density plots (a, b) and (c, d) respectively correspond to simulations carried with vorpal and impact-t. The beam energy is 
£ = 1 GeV. The DLW and bunch parameters are the ones appearing in Tab. [Til 



x^mVO ) ■ 

(25) 

In the latter equation the substitution cosh 4-> sinh 
should be made when considering dipole modes. 

The simulated electromagnetic field components are 
compared with vorpal simulations in Fig.[S]for the same 
case as presented in Fig. [51 In the Fig. [51 the impact- 
T simulation are performed with and without activating 
the space-charge algorithm. When accounting for space- 
charge forces, iMPACT-T is in very good agreement with 

VORPAL. 

The convolution over the z-summation in Eq. [531 
can be be performed with a numerical Fast Fourier 
Transformation (FFT). Therefore the total number 
of operations needed to evaluate the wakefields is 
oc NmodesN^NyNz log Nz. The total number of modes 
is twice the product between the modes allowed for each 



of the transverse wave numbers: Nmodes = '^Nk^Nk^- 
The factor of 2 comes from the inclusion of the dipole 
modes along with the vertically symmetric monopole 
modes. For the data generated in Fig. |8l the impact-t 
simulations are more than two orders of magnitude 
faster than the vorpal ones. 

The altered version of impact-t was used to explore 
a possible DLW experiment at the ASTA facility as a 1- 
GeV electron beam is injected in a DLW. The beam and 
structure have the parameters displayed in Tab. |TT] . For 
these simulations, the electron beam is an idealized cold 
beam without energy spread nor divergence. Figure [9] 
compares the longitudinal phase spaces simulated with 
VORPAL (top row) and impact-t (bottom row) at two 
axial locations. The agreement on the phase space struc- 
ture developing as this non-optimized beam propagates 
in the DLW is excellent. The mean and rms energies 
evolution for the witness and drive bunches reported in 
Fig. 1101 confirms the quantitative agreement between the 
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FDTD and semi-analytical models. Both models give 
mean and rms energies in agreement within ^ 5 %. 



1060 
1040 

> 

S 1020 



■ i 



Witness Beam 



il 



P 1000^" , , 

c 

^ 980 



960 



Drive Beam 



"^T) 2 4 6 8 10 12 14 16 18 20 22 

z (cm) 

FIG. 10: Evolution of mean (symbols) and rms (bars) energy 
as the drive (top traces) and witness (bottom traces) bunches 
propagates along the DLW structure. The simulations are 
carried with vorpal (red traces and symbols) and the modi- 
fied version of impact-t (blue traces and symbols). The beam 
energy is £ = 1 GeV. The DLW and bunch parameters are the 
ones gathered in Tab. [Ill The ordinate z = cm corresponds 
to the DLW entrance. 



VIII. CONCLUSIONS 

In this paper we presented a three-dimensional model 
to evaluate wakefields in slab- symmetric dielectric-lined 



waveguides. The only limitation of this model stems from 
the assumed charge distribution symmetry with respect 
to the vertical axis [p{x) — p{~'x)]. The model was 
successfully validated against three-dimensional FDTD 
PIC simulations performed with vorpal and was im- 
plemented in the popular beam dynamics tracking pro- 
gram iMPACT-T. The added capability to impact-t en- 
ables start-to-end simulation of linear accelerators based 
on DLW accelerating structures. Furthermore, because 
iMPACT-T includes a space charge algorithm, the up- 
graded version provides a valuable tool for investigating 
the performances of DLW acceleration when the dynam- 
ics of either, or both, of the drive and witness bunches is 
significantly impacted by space charge effects. The ob- 
served good agreement between the developed algorithm 
and simulations performed with the vorpal FDTD PIC 
program demonstrates that our model strikes an appro- 
priate balance between, efficiency, accuracy, and simplic- 
ity- 
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